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Abstract 

A model recently proposed for the calculation of air- fluorescence yield excited by electrons is revisited. Improved 
energy distributions of secondary electrons and a more realistic Monte Carlo simulation including some additional 
processes have allowed us to obtain more accurate results. The model is used to study in detail the relationship 
between fluorescence intensity and deposited energy in a wide range of primary energy (keVs - GeVs). In addition, 
predictions on the absolute value of the fluorescence efficiency in the absence of collisional quenching will be presented 
and compared with available experimental data. 
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1. Introduction 

The accuracy in the energy determination of 
ultra-high energy cosmic rays using the fluorescence 
technique is presently limited by the uncertainty 
in the air-fluorescence yield. In order to improve 
the accuracy of this fundamental parameter, sev- 
eral laboratory measurements are being performed 
|1|2|3|4|5|6|7|8| . On the other hand, progress on the 
theoretical understanding of the various processes 
leading to the air-fluorescence emission is being 
carried out. At high pressure (> 1 Torr) and high 
electron energy (> 1 keV) most of the fluorescence 
light is generated by low energy secondary elec- 
trons arising from N2 ionization [9J . A simple model 
which accounts for the contribution of secondary 
electrons to the fluorescence yield of the N2 2P sys- 
tem (C 3 II M — > B 3 II^) has been recently published 
[9]. Using the above model, the total fluorescence 
yield has been computed for the first time at high 
electron energies in [10] . 

Since fluorescence is mainly generated by sec- 



ondary electrons, their energy distribution is a key 
parameter. Previous results |9|10| were obtained us- 
ing the analytical energy distribution given by Opal 
et al. [11] which has been experimentally checked at 
low energy (below a few keVs). Unfortunately, the 
extrapolation of this formula to higher energies is 
not suitable to account for the energy distribution 
of delta rays properly. In the present paper, a new 
analytical approach valid in the whole energy inter- 
val (eV - GeV) is used. Preliminary updated results 
of our model using this energy spectrum together 
with improved molecular parameters and a more 
realistic Monte Carlo simulation are presented in 
this paper. 

Proportionality between fluorescence intensity 
and the energy deposited by the primary electron 
is usually assumed. Experimental results seem to 
support this assumption, at least, in small energy 
intervals. Recently the MACFLY collaboration has 
extended this experimental test to a large energy 
range, 1.5 MeV - 50 GeV, finding proportionality 
within their experimental errors of about 15% [3]. 
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However, from a theoretical point of view this pro- 
portionality is not obvious. The excitation of the 
fluorescence emission is strongly peaked at low en- 
ergies, in particular for the 2P system of nitrogen 
and therefore the ratio of fluorescence emission and 
deposited energy is strongly dependent on the spec- 
trum of low energy secondaries. In principle, this 
spectrum depends on both the primary energy and 
the distance from the primary interaction, and thus, 
this proportionality has to be demonstrated with a 
detailed analysis. 

2. The model 

The model described in |9)10| shows that the total 
number of photons of the v — v' band generated by 
an electron per unit path length in a thin target can 
be expressed by 

e vv ,(P) = N 1 + p /P , {<r V v>(E) + a vv ,(E,P)* ion (E)}. (1) 

N is the number of nitrogen molecules per unit 
volume. Collisional de-excitation (quenching) is 
taken into account by the well-known Stern- Volmer 
factor which depends on the characteristic pressure 
P' v . As shown in this expression, the fluorescence 
intensity is a sum of two contributions: a) direct 
excitation of v — v' band by the primary electron, 
given by the optical cross section a vv > and b) con- 
tributions from secondary electrons arising from 
molecular ionizations, proportional to the ioniza- 
tion cross section. 

The optical cross section is defined by 

Aw 1 / n \ 
a vv f = a v — — , (2) 

where a v is the cross section for electron excitation 
(ionization) of molecular nitrogen in the ground 
state to the upper v level of an excited electronic 
state of N 2 (Nj), and A vv > are the v — v' radia- 
tive transition probabilities (Einstein coefficients). 
Therefore the optical cross sections cr 2 p( _o), 
0"iiV(o-o) use d below are associated to the emis- 
sion of a 337 nm and 391 nm photon respectively 
induced by electron collision on a N 2 molecule. 

For high pressures (> 1 Torr) and energies 
(> 1 keV) the contribution from secondary electrons 
is dominant for generation of 2P fluorescence and 
relevant for the IN bands. The parameter a vv > in 
(PQ) represents the average number of v — v' photons 
generated inside the observation volume as a result 
of a primary ionization in the absence of quenching. 
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Fig. 1. Cross sections (atomic units) involved in the emission 
of air fluorescence as a function of energy. Cross section 
(Jexc accounts for all excitation processes without emission 
of a secondary electron. Cross section cr^ on accounts for all 
processes leading to the emission of a secondary electron. 
The cross section for the ejection of a K shell electron gk 
is shown for comparison. The optical cross sections for the 
excitation of the IN (0-0) and 2P (0-0) bands are also shown. 

The evaluation of a vv > has to be carried out by a 
Monte Carlo simulation which takes into account the 
various generations of secondary electrons produced 
after a primary ionization. The a vv > parameter turns 
out to be a function of primary energy E : pressure 
P and the geometrical features of the observation 
region. 

Next we will describe the molecular parameters 
used in our simulation. In the first place, figure 1 
shows the nitrogen cross section versus energy for 
various e~ -molecule interaction processes of inter- 
est in our problem. The cross section for molecu- 
lar excitation a exc accounts for all excitation pro- 
cesses without emission of a secondary electron. On 
the other hand, the ionization cross section <Ji on 
accounts for all processes leading to the emission 
of a secondary electron and can be described at 
large energy by the Born-Bethe function. At high 
energy only optically allowed transitions, also well 
described by a Born-Bethe function, contribute and 
therefore cr exc is proportional to cr ion . The cross sec- 
tion for the ejection of a K shell electron <jk is shown 
for comparison. Figure 1 also shows the optical cross 
sections (j2j) for the excitation of the 2P (0-0) and IN 
(0-0) bands of N 2 and N^ respectively. At low en- 
ergy, experimental values have been used (see [10] for 
references) and at high energies ctin(o-o) nas been 
extrapolated using the Born-Bethe function. For all 
cross sections the density correction at high energies 
|10|12|13] has been applied. 
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Fig. 2. Energy spectrum of secondary electrons emitted in a 
primary ionization induced by an electron of energy 1 MeV. 
The Opal formula is not able to reproduce the energy spec- 
trum of delta rays (Moller scattering). The extended Opal 
formula used in this work fits properly the Moller spectrum. 
The values of the EEDL database are also shown for com- 
parison. In ordinates the absolute value of the differential 
cross section in atomic units per eV is represented. 

Another ingredient of the model is the energy lost 
by the electrons in each collision process. For molec- 
ular excitation (without electron ejection) an av- 
erage value of < E exc > = 8.5 eV has been esti- 
mated. In an ionization process the lost kinetic en- 
ergy (primary and secondary electrons) equals the 
N2 ionization potential, I = 15.5 eV, plus the excita- 
tion energy of the ion, for which an average value of 
< E l °™ c > = 1.3 eV has been estimated at low en- 
ergies from available data. At large energy K-shell 
ionization is also taken into account with a relative 
probability given by the corresponding cross section. 
In this case the kinetic energy is reduced by 410 eV 
(K-shell binding energy). As a result < E l ° x n c > in- 
creases up to 5 eV. 

In the present calculations the tracks of secondary 
electrons have been simulated using available angu- 
lar cross section for the various individual interac- 
tion processes, including elastic e~-molecule colli- 
sions. Other processes like generation of X ray from 
molecular de-excitation and bremsstrahlung by high 
energy electrons have also been included in our sim- 
ulation although the final results are nearly inde- 
pendent on these processes. 

As already mentioned, the total fluorescence 
emission is very sensitive to the energy spec- 
trum of secondary electrons. At low energy the 
analytical expression proposed by Opal et al. 
fits accurately experimental results. At high pri- 
mary energies, for which no experimental data 
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Fig. 3. Same as figure 2 for 1 GeV primary energy. 

are published, precise calculations are described 
by Seltzer [12 13J. Numerical results for any atom 

HI 



are available in EEDL databases |13)14j . Instead 
of implementing these databases in our simula- 
tion, we looked for an analytical expression for 
the energy spectrum in the whole energy range 
for both primary and secondary electrons. The 
desired function must fulfill several requirements: 
a) it has to behave as the Opal formula at low 
primary energy; b) at high E the function has to ac- 
count for the well known energy spectrum of delta 
rays, c) the integral of the differential cross section 
has to give the total ionization cross section and 
d) the assumed function must be consistent with 
the collision energy loss given by the Bethe-Bloch 
theory for electrons. 

In this work a new analytical formula has been 
used which properly accounts for the above four con- 
ditions. Details on this analytical approach will be 
given in a paper in preparation. In figures 2 and 
3 the energy spectrum of secondary electrons pre- 
dicted by this "extended" Opal formula used in the 
work is compared with other available data for pri- 
mary electrons of 1 MeV and 1 GeV respectively. 

Fulfilment of requirement d) can be easily checked 
as follows. The energy loss per unit electron path 
length in air, assuming N air molecules per unit vol- 
ume, is closely related with the mean value of the 
energy distribution of secondary electrons by the fol- 
lowing relation 



dE 

dx 



N atr {< E° dep > + <E s >}a ion (E) , 



(3) 



where < E s > is the mean value of the energy dis- 
tribution of secondary electrons at a given E value 
and 
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Fig. 4. Energy loss of electrons per unit path length versus 
primary energy. Continuous line represents the calculated 
values from equation (|3} using the average energy of the 
extended Opal formula. Black dots are the values from the 
Bethe-Bloch formula. In both cases the medium is air at 
atmospheric pressure. 



< Kep > = < E e 



(4) 



represents the average value of the energy deposited 
in the medium by the primary electron per primary 
ionization process. 

In Figure 4 the energy loss calculated from ex- 
pression (|3]) using the < E s > value from our energy 
distribution of secondary electrons is represented, 
showing full agreement with the Bethe-Bloch for- 
mula [T5] . 

3. Results 

In this section we will present results on the var- 
ious parameters of interest for the interpretation of 
the air fluorescence measurements. In our simula- 
tion the observation region has been assumed to be 
a sphere of radius R with center in the interaction 
point of the primary electron. Assuming a fixed tem- 
perature, both fluorescence intensity and deposited 
energy turn out to be basically dependent on the 
product PR, that is, on the number of interactions 
inside the observation region. In fact, due to the den- 
sity correction at very large energies, for a given PR 
value, the above parameters are slightly dependent 
on pressure. However, the corresponding deviation 
is not significant inside the range 100 - 760 Torr. 

The energy deposited by the primary electron per 
unit path length can be calculated from the average 
energy deposited by secondary electrons < Ed ep > 
as 
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Fig. 5. Energy deposited by a primary electron per unit 
path length versus primary energy for several values of PR 
(continuous lines). For very high PR value deposited energy 
equals the energy loss of the primary electron predicted by 
the Bethe-Bloch theory (broken line). 



dx 



N a i r {< > + < Edep >}cion(E) . 



(5) 



Note that this equation is the same as (j3]) replac- 
ing < E s > by < E dep >. The value of < E dep > 
has been obtained from our Monte Carlo simulation. 
Figure 5 shows the energy deposited per unit path 
length versus primary energy for several PR values. 
We have checked that, as expected, the deposited 
energy for an unlimited medium, PR — > oo, equals 
the energy loss predicted by the Bethe-Bloch theory. 

The fluorescence efficiency & vv >, defined as the 
ratio between the number of fluorescence photons 
emitted in a given v — v' molecular band and the 
deposited energy is a fundamental parameter. Com- 
bining equations $1]) and ([5]), the following useful 
relation can be obtained: 



vv' 



N 

-^V air 



< E l P > 



< Edep > 



(6) 



where $® v , is the fluorescence efficiency in the ab- 
sence of collisional quenching and thus related with 
®w' by 



(7) 



■P/Pl 



The values of a vv > and $>® v , have been computed 
for the 2P (0-0) and IN (0-0) bands, for which ac- 
curate experimental results of the corresponding 
cross sections are available. Figures 6 and 7 show 
$° against E for several values of PR while in fig- 
ures 8 and 9 the dependence of <l> versus PR is 
represented for various energies. 

These figures show several interesting features. 
In the first place, for high energy (beyond the keV 
range), and for observation regions larger than 
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Fig. 6. Fluorescence efficiency for the 2P (0-0) band in the 
absence of collisional quenching against primary energy for 
several values of PR (continuous lines). Circles represent the 
result for a very large medium. 
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Fig. 7. Same as figure 6 for the IN (0-0) transition. 

about 1 Torrxcm, <£° is nearly constant. Never- 
theless a smooth dependence on primary energy is 
found. The value of 3>° decreases with primary en- 
ergy about 10% in the range 1 keV - 1 MeV and 4% 
in the interval 1 MeV - 20 GeV for the 337 nm band 
(figure 6). For the 391 band the corresponding de- 
creasing is about 6% for the interval 1 keV - 1 MeV 
and 1% for 1 MeV - 20 GeV (figure 7). On the other 
hand a smooth growing of <£° with PR, smaller than 
2% in the range 10 - 1000 Torrxcm, is also found 
for energies larger than 1 MeV (figures 8 and 9). At 
lower energy and/or region size the fluorescence is 
clearly not proportional to the deposited energy. 

The real number of photons detected per unit of 
deposited energy <£> can be calculated from equation 
Q using the corresponding P' v value. Unfortunately 
there are significant disagreements in the available 
experimental data of collisional quenching in air. 
Therefore the comparison of the absolute value of 
the efficiency with experimental data is not straight- 
forward. On the other hand, some experiments only 
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Fig. 8. Fluorescence efficiency for the 2P (0-0) band in the 
absence of collisional quenching against the size of the ob- 
servation region PR for several primary energies. 
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Fig. 9. Same as figure 8 for the IN (0-0) transition. 

provide a measure of the total fluorescence in a wide 
spectral range including most of the molecular tran- 
sitions. In addition, very often the fluorescence yield 
is given in units of photons per meter. 

The fluorescence efficiency at zero pressure for a 
reference molecular band (e.g. the 337 nm band) can 
be inferred from the number of photons per unit elec- 
tron path length e\ emitted in a given wavelength 
interval at a given pressure by 



^337 



h37 



dE dep /dx I\ 



P' 

^337 



(8) 



where the ratio between the fluorescence intensity 
for the 337 nm band ^337 and that for a broad 
wavelength interval I\ can be calculated from the 
Franck- Condon factors, Einstein coefficients and 
P' v values, as explained in detailed in ref. [TP] , For 
the comparisons shown below the energy deposited 
per unit path length dEdep/dx has been obtained 
from our MC simulations (figure 5) assuming for all 
experiments a typical PR value of 1600 Torrxcm. 
Notice that this parameter depends on PR very 
smoothly (around 6% variation in the interval 300 - 



5 



Table 1 

Parameters for the calculation of ^§37 experimental data at atmospheric pressure. Experiment reference, electron energy, 
spectral bandwidths and measured number of air-fluorescence photons per meter are displayed in columns one to four. Fifth 
column shows the deposited energy per meter obtained from our simulations. The sixth one gives the total fluorescence 
efficiency at atmospheric pressure in the spectral interval. The ratio of the 337 nm band intensity to total intensity is shown 
in next two columns (7 and 8). Ratios have been calculated assuming the set of P' v values given either by references [16] (P' p 
column) or [5 [P' A column). Next two columns show inferred values of ^§37 assuming the above mentioned P' v values. Column 
eleven shows the ^§37 result reported by the experiment. Last column shows the preliminary theoretical results presented in 
this work. Bolded figures are those inferred from our calculations. 



Experiment 
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[nm] 


£\ 


dE dep /dx 
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r [MeV]- ] 
reported 


this work 


Nagano et al. 
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17.6 


0.27 


0.26 


380 
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t Value reported by FLASH. 



3000 Torr x cm) . For the following discussion we will 
use P' v measurements very recently published by the 
AIRFLY collaboration [5J and those obtained from 
deactivation rate constants N-N and N-0 measured 
by Pancheshnyi et al. [16J (see also [TO]). 

We have compared our theoretical §% 37 value 
with some available experimental data. For instance 
Nagano et al. PQ, working at an average energy 
of about 0.85 MeV, reports a value of 272 MeV" 1 
significantly smaller than our predictions. However 
from their measurement of 1.021 photons (337 nm) 
per meter at atmospheric pressure, results from (|8]) 
of 440 or 370 MeV -1 are found by using P' v values 
of [16] or [5] respectively, the last one in very good 
agreement with our predicted value of 365 MeV -1 
(see table 1). This discrepancy can be explained by 
comparing the P 337 value of 14.4 Torr reported by 
PQ with those of 9.8 and 11.9 Torr from 16] and 
[5] respectively. This simple exercise shows that the 
result of the comparison is strongly dependent on 
the chosen P 337 value. 

Another interesting example is the comparison 
with the measurement recently reported by the 
FLASH collaboration at 28 GeV [3]. In this case 
we use their 5.06 photons/m for the whole fluores- 
cence spectrum in the wavelength range 300 - 420 
nm. The contribution to the total intensity I\ of 
the Gaydon-Herman bands measured by AIRFLY 
[5] has been included. The inferred $337 results are 
453 and 353 MeV -1 using P 3 / 37 values of [16] or [5] 
respectively (see table 1). Again the second one is 
in very good agreement with our predictions. Note 
that, as shown in the table, the factor I337/IX is 
only weakly dependent on the P' v set used. 

Finally the $g 37 result reported by MACFLY [4] 



of 174 MeV -1 , constant in the range 1.5 MeV - 50 
GeV, is smaller than our predictions. Notice that 
the value P 337 = 19.4 Torr of [4] is larger than those 
of [16] or [5]. In fact, following the above described 
procedure a $337 value of 380 or 290 MeV -1 , using 
[T6] or [5] respectively, would be inferred from the 
MACFLY result of $ A = 17.6 photons/MeV at at- 
mospheric pressure in the spectral interval 290 - 440 
nm. 

4. Conclusions 

An improved model for the air-fluorescence emis- 
sion has been used to analyze in detail the relation- 
ship between fluorescence intensity and deposited 
energy. For primary energies larger than a few keV 
and assuming observation regions over 10 Torr x cm 
nitrogen fluorescence is nearly proportional to de- 
posited energy. A smooth decrease of of about 
10% in the range 1 keV - 1 MeV and 4% in the in- 
terval 1 MeV - 10 GeV for the 2P system is found. 
Taking into account that the contribution of the en- 
ergy released by an extensive air shower by electrons 
with energy smaller than 1 MeV is only of about 
22% [17 , the above smooth dependence of efficiency 
on primary energy has no impact on the calibration 
of fluorescence telescopes. Nevertheless very precise 
laboratory measurements (better than 10%) of the 
fluorescence yield in the range keVs - GeVs could 
observe this small effect. 

The absolute value of the fluorescence efficiency 
for the 337 nm band computed in this work has 
been compared with some experimental results. The 
comparison was carried out by reducing experimen- 
tal values of number of photons per meter at atmo- 
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spheric pressure to ^ 37 using quenching data from 
other authors. In general we found a good agree- 
ment, in particular when using quenching parame- 
ters from AIRFLY. 

Our improved energy distribution of secondaries 
(consistent with the energy loss predicted by the 
Bethe-Bloch formula) has lead us to demonstrate 
quantitatively that fluorescence intensity is basi- 
cally proportional to deposited energy. On the other 
hand, uncertainties in the absolute values of theo- 
retically predicted &® v , quantities are larger because 
these are limited by the accuracy of the assumed 
molecular parameters involved in the calculations. 
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